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We study the strong-coupling superconductivity near ferromagnetic quantum critical points, mainly fo- 
cusing on the upper critical fields H C 2- Based on our simple model calculations, we discuss experimentally 
observed unusual behaviors of H C 2 in a recently discovered ferromagnetic superconductor UCoGe. Espe- 
cially, the large anisotropy between H C 2 || a-axis and H C 2 || c-axis, and the strong-coupling behaviors in 
H^2 are investigated. We also examine effects of non-analytic corrections in the spin susceptibility on the 
superconductivity, which can arise from effective long range interactions due to particle-hole excitations. 

KEYWORDS: superconductivity, upper critical field, ferromagnetism, quantum criticality, non-analytic 
corrections 



1. Introduction 

Recently, a new ferromagnetic superconductor UCoGe was 
discovered and has been attracting much interest. 1 " 8 ' It is con- 
sidered to be located near a ferromagnetic quantum critical 
point (QCP), and the related phenomena are extensively stud- 
ied. In this compound, the superconductivity coexists with 
ferromagnetism in some pressure range, and the pairing state 
is expected to be spin triplet with equal-spin pairing along the 
easy axis of the magnetization, c-axis. The superconductivity 
extends to the paramagnetic phase and the superconducting 
transition temperature changes smoothly over the two phases, 
which suggests that the superconductivity is mediated by the 
Ising ferromagnetic spin fluctuations. 9 " 11 ' 

Intensive experimental studies reveal interesting behaviors 
of the upper critical field H c2 which would be related to 
quantum criticality. 2 ' 4 ' 5 ' Naively, it may be expected that the 
anisotropy in the upper critical field H c2 in a superconduc- 
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tor with equal-spin pairing obeys H" 2 > H££, because 
there is no Pauli depairing effect for c-axis. In the coexis- 
tence region of the ferromagnetism and the superconductiv- 
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ity in UCoGe, however, the observed H c2 is H^f 3> H t 



which is completely different from the naive expectation. Re- 
cently, it was pointed out that, for the superconductivity coex- 
isting with the ferromagnetism which has a sufficiently large 
magnetization, the Pauli depairing effect for a field perpen- 
dicular to the magnetization does not destroy the supercon- 
ductivity and the orbital depairing effect would be essential 
there. 12 ' Although this explains a very important part of the 
interesting behaviors of H c2 in UCoGe, there remain unre- 
solved problems which would be beyond the weak-coupling 
theory. In addition to the unusual anisotropy, the behaviors 
of H c2 || a, b are surprising. 4,5 ' H'^'s are huge exceeding 
10 (T) while the transition temperature T sc at zero field is 
0.8 (K) and have upward curvatures, which would be con- 
sidered as a characteristic property of the strong-coupling su- 
perconductivity. The huge H c2 with upward curvatures are 
also observed in noncentrosymmetric heavy fermion super- 
conductors Ce(Rh/Ir)Si3, 13 ' 14 ' which can be understood as an 
interplay of lack of inversion symmetry and quantum critical- 
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ity. 15 ' For H^ 2 in UCoGe, H c2 shows an S-shaped curve at 
high fields and the similarity to that in URhGe was pointed 
out. 5 ' These behaviors would be related to the ferromagnetic 
quantum criticality in UCoGe. 

Apart from experiments, there have been new theoretical 
progress on ferromagnetic quantum criticality. 16 " 19 ' It is ar- 
gued that the criticality can be affected by effective long range 
interactions due to particle-hole excitations which lead to the 
so-called non-analytic corrections. If the superconductivity is 
mediated by ferromagnetic spin fluctuations, it could be also 
affected by them. 

In this study, we examine the strong-coupling superconduc- 
tivity near ferromagnetic QCPs having UCoGe in mind. We 
restrict our calculations to the paramagnetic states and ne- 
glect the Pauli depairing effect. Under these conditions, we 



show that h\ 2 ^> H^ 2 can hold in some parameter range. We 
perform the same calculations for two cases, the case of the 
point node symmetry and the line node symmetry, and con- 
clude that the point node symmetry is a promising candidate 
for the superconductivity realized in UCoGe. We also study 
the effects of non-analytic corrections in the spin fluctuations 
on the superconductivity. 

2. Calculations and Results 

In this study, we use a very simple model to investigate 
the superconductivity near the ferromagnetic QCP. 9 " 11 ' The 
action is 
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Sq = 

H, and the 



Efc4+ g ( CT V2)cfe, e k 
filling is fixed at n — 0.15 for which the Fermi surface 
is almost a sphere. We have taken the energy unit t = 1 
and the length unit a = 1. The spin fluctuations are of 
Ising-type according to the NMR experiments 8 ' and the re- 
sulting pressure-temperature phase diagram can be consis- 
tent with the experiments. 9 " 11 ' The criticality is characterized 
by S = 5q(T + 9) + Ch(nBH z ) 2 , where we have just as- 
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sumed the mean-field-like temperature dependence. Here, we 
fix So = 4 and v = 4. The important point is that the spin fluc- 
tuations depend on the applied magnetic field along the c-axis 
H z , as suggested by the resistivity measurements 5 ' and the 
NMR experiments. 8 ' The non-analytic correction is incorpo- 
rated with the coefficient c„ which can arise from the effective 
long range interactions due to particle-hole excitations. 1 6-1 9 ' 
If we define the effective mass of the fluctuating modes as 
<5 e fi = S + min[c„g 2 \nq + q 2 ], S e s becomes smaller than S 
for c n > 0, which means that the spin fluctuations are en- 
hanced by the non-analytic correction. Here, we just put a re- 
mark on how actually the non-analytic correction can appear. 
In a fully SU(2) symmetric one-band model, all the scatter- 
ing processes with momentum transfer q ~ cancel out and 
cannot contribute to the spin susceptibility because of the loop 
cancellation theorem. 20 ' For such a case, the non-analytic cor- 
rection can arise from the scattering processes which include 
at least one scattering with q ~ 2fcp, where kp is the Fermi 
wavenumber. On the other hand, the system with Ising-like 
symmetry, there would be no strong reason by which q ~ 
scattering processes completely cancel out. This issue would 
be discussed in detail elsewhere. 

To study the superconductivity focusing on the orbital de- 
pairing effect, we solve the Eliashberg equation without the 
Pauli depairing effect. The linearized Eliashberg equation 
reads, 

A aa (k) = -^]Ty(fc,fcW)A^(fc'), (3) 

k' 

g(k) = i[(^ o |G(fc + n)G(-fc)|0o) 

+(<h\G(k)G(-k + U)\<f>o)] (4) 

where n = (0, -iV R - 2eA(R)) with V x A = (H, 0, 0) 
or (0, 0, H), and R is the center of mass coordinate of the 
Cooper pair. </>o is the lowest level Landau function and G 
is the Green's function with the selfenergy evaluated at the 
lowest order in fl 2 Xo ,S(fc) = (T/3N)J2 q g 2 X (q)G (k + q), 
where Go is the non-interacting Green's function. The pairing 
interaction is also calculated at the lowest order, V(k, k') = 
-(l/6)3 2 xO - k ') + (l/6).9 2 x(^ + k'). The calculations 
are similar to those in the previous study. 15 ' We fix t = 50 
(K) and a = 4.0 (A) in this study, which corresponds to the 
effective mass of the cyclotron motion of the orbital depairing 
m cff = h 2 /(ta 2 ) ~ 440 x (bare electron mass). 
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Fig. 1. (Left panel) Transition temperature T sc as a function of 9 for c n = 
0. (Right panel) T sc as a function of c n for 6 = 0. 

First, we study the effects of the non-analytic correction in 
the spin fluctuations on the superconductivity at zero field. 
At zero field, the transition temperature for the gap func- 
tion d = (sin/c a ,sinfc{,,0) with the point nodes and the 



gap function d = (sinfc c ,0,0) with the horizontal node are 
degenerate. In this study, the coupling constant is fixed as 
g 2 Xo = 112. Qt which is large so that we have rather high 
transition temperatures within our model. For this value of 
g 2 Xo, the mass enhancement factor due to the ferromagnetic 
spin fluctuations is z ~ 0.8. The left panel of Fig. 1 shows the 
^-dependence of the transition temperature T sc which corre- 
sponds to the pressure-temperature phase diagram of the para- 
magnetic side in the experiment and 6 = corresponds to the 
QCP. The transition temperature T sc is highest at 9 = and 
is decreased as 9 increases, which is consistent with the ex- 
perimental phase diagram of UCoGe. 3 ' 4 ' In the right panel of 
Fig. 1, we show T sc as a function of the coefficient of the 
non-analytic correction c n . We see that T sc is monotonically 
enhanced as c„ increases. This is because the effective mass 
of the criticality S c g becomes small with finite c„ > 0. This 
shows that the enhanced criticality by the non-analytic cor- 
rection favors the superconductivity in the case of Ising-type 
spin fluctuations. 

Next, we move to the effects of a magnetic field on the su- 
perconductivity and fix 9 = 0. In the presence of the mag- 
netic field, we use the same model, eq.(2), although there 
might be other non-analytic corrections related to the mag- 
netic field. 19,21 ' Even if we neglect Pauli depairing effects, the 
degeneracy of the symmetries of the gap functions is lifted by 
H and the anisotropy of H c2 depends on the positions of the 
gap nodes. 22 ' This is simply explained as follows. The effec- 
tive velocity for the cyclotron motion of the orbital depairing 
effect can be of the form of Vk = (p(k)vk for the basis func- 
tion <p corresponding to the superconducting gap symmetry. 
When the magnetic field is parallel to the c-axis, the cyclotron 
motion can enjoy the nodes for the case with horizontal line 
nodes where Vk — 0, resulting in a large orbital limiting field. 

The cyclotron motion cannot do so for H _L c, which leads 

He 

to the anisotropy of the orbital limiting field H" 2 > H£f. On 
the other hand, under H \\ c, Vk cannot be zero for the point 
node gap function except for the poles of the Fermi surface, 
k a = kfy = 0. Therefore, for point nodes at the poles of the 
Fermi surface, the order is turned over, < H^ c . 






H||a ■ 




H]]c — 







0.25 0.5 



Fig. 2. Temperature v.s. upper critical fields H c 2 for c„ = (left panel) 
and c n = 0.4 (right panel) for d = (sin k a , sin k b , 0). The solid curves 
are for H c 2 \\ a-axis, and the dotted curves are for H c 2 \\ c-axis with 
Ch = 0, Ch = 1.0, Ch = 10.0 from the top to the bottom, respectively. 

In Fig. 2, we show the upper critical field for several val- 
ues of c/j for the point node case for which d-vector is d = 
(sin k a , sin kb, 0). H c2 is normalized by the transition temper- 
ature at zero field T sc0 which is a measure of the Pauli lim- 
iting field for the usual s-wave superconductivity. The solid 
curve is for H \\ a and dotted curves are for H \\ c. We see 
that, even if Ch = 0, > holds because of the point 
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nodes. When the Ising spin fluctuations are suppressed by the 
magnetic field H \\ c as suggested by the experiments, 5 ' 8) 

is reduced and it has only moderate temperature depen- 
dence for sufficiently large Ch- If there exists strong suppres- 
sion of the spin fluctuations, H c2 can have strong anisotropy 

^> H^- Note that, for the Ising anisotropic case, the spin 
fluctuations are robust against in-plane fields, resulting in that 

remains large with the strong-coupling behaviors. The 
non-analytic correction can make the result more drastic and 
the calculated anisotropy is quite strong as seen in the right 
panel of Fig. 2. 

The strong enhancement of the calculated is under- 
stood as a result of an increasing pairing interaction and a de- 
creasing depairing effect of the spin fluctuations at low tem- 
peratures. The physical origin is the same as that of the huge 

lie 

H" 2 in noncentrosymmetric heavy fermion superconductors 
Ce(Rh/Ir)Si3. 13 ~ 15) Therefore, as pointed out in the previous 
work, 15 ' it can be said that the strong enhancement in the or- 
bital limiting field near the quantum criticality has universal 
nature. In this study, however, we have assumed that the su- 
perconductivity is orbital limited, and the suppression of the 
Pauli depairing effect 12 ' would work well in the coexistence 
region of the ferromagnetism and the superconductivity in 
UCoGe. In the coexistence region, the mass of the criticality 
S is replaced with the magnetization which has weak temper- 
ature dependence. Whether such a mass can actually lead to 
huge without the assumption of the orbital-limited su- 
perconductivity is an open issue which should be addressed 
in the near future. 

On the other hand, for the case of the gap function d = 
(sin k c , 0, 0) with a horizontal line node, the relation < 
H^2 holds when the suppression of the spin fluctuations is 
not so strong as shown in Fig. 3. For H c2 to be > H^, 
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Fig. 3. Temperature v.s. upper critical fields H c i for c n = (left panel) 
and c„ = 0.4 (right panel) for d = (sinfc c ,0, 0). The solid curves are 
for H c 2 || a-axis, and the dotted curves are for H c 2 || c-axis with c h = 
0, Ch = 1.0, = 10.0 from the top to the bottom, respectively 

a quite large c/, is needed, which cannot be usually realized. 

From the above results, we conclude that a gap function 
with point nodes is more promising than that with horizon- 
tal line nodes as a candidate for UCoGe in the coexistence 
region. In UCoGe which has a orthorhombic structure, possi- 
ble gap functions for the coexistent states are classified by the 
magnetic point group £>2(Cf ). 23) Among the possible sym- 
metries, the most promising candidate is the Ai state or the 
A2 state for H C 2 to be consistent with the experiments. 2 ' 4 ' 5) 

In summary, we have discussed the superconductivity near 
the ferromagnetic QCP within a simple model for the Ising 
spin fluctuations. The non-analytic correction can enhance 



the superconducting transition temperature. Under a magnetic 
field, the positions of the nodes strongly affect the orbital 
cyclotron motion, leading to the anisotropy in H c2 . Taking 
into account the suppression of the Ising spin fluctuations by 
H || c, we have shown that the anisotropy can be <C 
for the point node case, and has characteristic nature 
for strong-coupling superconductivity. On the other hand, it 
is rather hard to reproduce such behaviors for the horizontal 
line node case. These results suggest that the superconductiv- 
ity realized in UCoGe is in the Ai state or the A2 state, both 
of which have no horizontal line nodes. 
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